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Abstract 

Laser-plasma interferograms are currently analyzed by extracting 
the phase-shift map with FFT techniques (K.A.Nugent, Applied Op- 
tics 18, 3101 (1985)). This methodology works well when interfer- 
ograms are only marginally affected by noise and reduction of fringe 
visibility, but it can fail in producing accurate phase-shifts maps when 
dealing with low-quality images. 

In this paper we will present a novel procedure for the phase-shift 
map computation which makes an extensive use of the Ridge Extrac- 
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tion in the Continuous Wavelet Transform (CWT) framework. The 
CWT tool is flexible because of the wide adaptability of the analyz- 
ing basis and it can be very accurate because of the intrinsic noise 
reduction in the Ridge Extraction. 

A comparative analysis of the accuracy performances of the new 
tool and the FFT-based one shows that the CWT-based tool phase 
maps are considerably less noisy and it can better resolve local inho- 
mogeneties. 

1 Introduction 

Interferometric techniques are widely used to characterize the optical prop- 
erties of a variety of media. An important class of applications concerns the 
investigation of the density distribution of plasmas produced by high inten- 
sity laser-matter interactions. In recent years various interferometer schemes 
have been developed and successfully applied to the characterization of the 
wide range of plasma condition which can be achieved in laser-plasma ex- 
periments, from the long-scalelength underdense plasma generated by laser 
explosion of a thin foil target to the steep, denser plasma generated by short 
pulse interaction with a solid target. All these schemes make use of a so- 
called probe beam which consists of a laser pulse which probes the plasma at 
a given time IQ 0. The fringe pattern must be then analysed to obtain the 
two-dimensional phase-shift which contains the physical information on the 
plasma. Then, provided that appropriate symmetry conditions are satisfied, 
inversion techniques can be applied to the phase shift map to obtain the 
density profile 0. 

The traditional way of reading a fringe pattern consisted in building a 
grid over the pattern and in measuring, for each position on the grid, the 
number of fringe jumps with respect to the unperturbed fringe structure. 
This procedure was very simple and was performed manually. However, the 
amount of information which could be extracted in this way was very limited 
due to the small number of grid points that can be employed. 

In 1982 a novel fringe analysis technique was proposed in which the 
phase extraction was carried out using a procedure based upon Fast Fourier 
transform (FFT). This technique allows the information carried by the fringe 
pattern to be decoupled by spatial variations of the background intensity as 
well as by variations in the fringe visibility, provided that the scalelength of 
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such perturbations is large compared to the fringe separation. A few years 
later this FFT technique was applied for the first time to laser-produced 
plasmas 0. More recently the technique was extensively applied by our 
group to the analysis of long-scalelength underdense laser-plasmas . The 
use of ultra-short probe pulses has considerably reduced the fringe-smearing 
effect due to the motion of the plasma during the probe pulse. This fact 
allowed to investigate short-lived phenomena in the propagation of ultra- 
short laser pulses through plasmas [§]. 

The extensive use of the FFT based technique carried out by our group 
has shown its effectiveness. In some circumstances however, reduction of 
fringe visibility, non uniform illumination, noise and the presence of local 
image defects make the FFT based technique unstable and the results are 
not fully satisfactory because of the presence of unphysical phase jumps. In 
addition small scale non-uniformities with low departure from the density 
background are unlikely to be detected. 

In this paper we show that Continuous Wavelet Transforms can also be 
applied to the analysis of interferograms resulting in a more flexible and 
reliable technique than the FFT based one. To our knowledge, this is the 
first time that such an approach is applied to fringe pattern analysis. 

In section 2 we will shortly introduce the Continuous Wavelet Transform 
(CWT) and its remarkable properties of good space-scale analyzer. 

In section 3 we introduce our I ACRE, " Interferogram Analysis by Con- 
tinuous wavelet trasform Ridge Extraction" tool for interferograms analysis 
and we compare its sensitivity to the FFT-based method. 

Section 4 is devoted to conclusions and comments. 

2 The Continuous Wavelet Transform analy- 
sis tool 

The Continuous Wavelet Transform is a tool to obtain a representation of 
signal s which is intermediate between the "real time" description s = s{t) 
and the "spectral" description s = s{uj), so that it is a very powerful tool 
to obtain a time-frequency description of a sequence of data. In this paper 
the words "time" and "space", or "frequency" and "spatial-frequency" will 
be used indifferently. 

The need of a time-frequency description of a sequence is much strong 



3 



when the signal represents a sum of frequency-modulated components (as for 
each section of an interferogram image, as shown below). This is because in 
a purely spectral analysis the frequency content of a modulated sinusoid is 
generally spread in a large region and no identification of the signal from its 
spectral amplitude is allowable. 

The obvious step that can be made to overcome the lack of time sensitivity 
is the introduction of a sequence of windows of a given width and centered 
at different times: for each window the FFT of the signal is computed and 
a partial time resolution is obtained. These techniques are called "Short- 
Time Fourier Transform" or "Gabor Transform" |^. The Gabor Transform is 
currently used in many context but is not considered by the signal-processing 
community a "full analysis tool". This is because the number of oscillations 
of each sinus in the window depends on the frequency and consequently the 
spectral and spatial resolutions should be optimized (by tuning the window 
length) only in a narrow band. 

From the early 80 's, with the introduction of the Wavelet Transform, a 
satisfying time- frequency analysis tool is available [|11[] |T2| . 

To introduce the Wavelet Transform, let us first define notations for the 
Fourier Transform. For a signal s e L^{TZ)r\L'^{Tl) the Fourier coefficients, 
that is the scalar product between the signal and the infinitely oscillating 
terms = exp(— i cut): 



form a complete basis of the space to which s belongs. 

Let us introduce a function \l/(t) called Mother wavelet. Now, instead of 
decomposing the signal s as a sum of the pure oscillating terms e^^ {Fourier 
Transform), we build a decomposition of s in terms of the base of all the 
translated (by parameter b) and scaled (by parameter a) \l''s. The base of 
the Continuous Wavelet Transform (CWT) is then a two-parameter family 
of functions 



The choice of the Mother Wavelet used to build the analyzing base is 
quite free and must be adapted to the actual information that should be 



Once the base has been built, one can compute the CWT coefficients as 




(1) 





extracted from the signal |IT] . 
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the scalar product of the signal and '^a,b- 



Ws{a,b) = <^a,6|^^> (3) 



J-oo a \ a J 



A particular choice of Mother Wavelet is the Morlet wavelet, and is largely 
used in studying signals with strong components of pure sinus or modulated 
sinusoids. The Morlet base has the form 

^{t) = exp{iuot) exp(-(t/T)') , (4) 

where the parameters ujo and r control the peak frequency and the width 
of the wave respectively. The product cuq x t controls the time and spectral 
resolution of the Wavelet decomposition: a large r corresponds to a long wave 
(high spectral resolution and low temporal resolution) while a small r pro- 
duces an "event based" analysis (low spectral resolution and high temporal 
resolution). 

We now face the problem of a numerical computation of the Wavelet 
coefficients map Ws{a,b). For a sequence of N samples sf, i = 1...N of s, 
the translation parameter b (which controls the central position of the wave 
envelope) can be sampled in a straightforward way: b ^ bf, i = 1...N. The 
scaling parameter a (which controls the characteristic scale of the wave) may 
be sampled as aj — 2"^/^" j — 1...M, (Natural or Log sampling), where 
Ny is the "number of voices per octave" parameter. Each aj is called "voice" 
and, in the case A^^, = 12, Log sampling exactly corresponds to the spectral 
sampling of musical tones in the "tempered scale" introduced by J.S. Bach. 

The Log sampling of CWT coefficients in the Morlet basis is very useful 
when the spectral content of the signal is the main information to be ex- 
tracted, because it provides a good compromise between spatial and spectral 
resolution. As the reader can easily check, the spectral resolution at each 
voice is proportional to the peak frequency of the voice (j^^ so that the 

relative spectral uncertainty ^ is constant along the a axes. 

The real part of CWT map shows an important feature of CWT with the 
Morlet base: Tl{Ws) is almost constant, apart from the thin band centered 
at the local signal frequency. Futhermore, the sequence 

R,{b) = n{W,){b,aR{b)) (5) 
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where, for each b*, a{i{b*) is the voice corresponding to a local maximum 
of the the line-out of the absolute value CWT map taken at 6 = 6*, well 
reproduces the input signal itself. The sequence (ore more generally the 
sequences when more complex signals are analyzed) (H) is called the Ridge of 
CWT map and represents the subset of CWT map where most of the " energy" 
is contained. Presently, the Ridge detection of the CWT map plays a rising 
role in signal processing [0[|14|; especially in the search of non-stationary 



signals with a very low signal-to- noise ratio (see and references therein). 
This is because the Ridge sub-map well captures the "true" input signal 
even in the presence of a quite strong noise. Futhermore, Ridge extraction 
in CWT maps of analytic signals represents a natural way to detect the local 
frequency evolution and, eventually, to easily recover phase information. 



3 The new CWT-based method. 

3.1 The FFT-based method for phase-shift estimation 

The extraction of phase-shift map, that is the computation for each pixel 
of an interferogram image of the phase-shift with respect to a unperturbed 
wave profile, is usually performed with the help of Fast Fourier Transform 
(FFT-based method). Consider for example the interferogram of Fig. ^ of 
a laser-exploded foil target plasma 0. Let its gray-level map be I{z,x) and 
for each Z build the sequence sz = I{z = Z, x) (that is a horizontal line-out 
of the figure). The background fringe pattern would give sequences sz very 
similar to pure oscillating terms plus noise and (possibly) a slowly varying 
background. If the departure of such a behaviour is identified as a local 
frequency modulation of the oscillating term, then the phase-shift 6(j){z, x) 
can be easily computed as the difference between the perturbed phase at each 
X position and the corresponding phase of the background sequence. Figure § 
shows a sequence sz for Z = 400 (the middle of the frame). The behaviour of 
Sz can be identified as a frequency-modulated oscillation with local frequency 
fl{x) increasing with x, plus noise and slowly rising background. In addition, 
the amplitude of oscillations sharply reduces for x ~ 700 (this phenomenon 
is known as "reduction of fringe visibility", see |0). 

The FFT-based phase-shift extraction uses FFT for both filtering the 
sequence from noise and background (with cuts in the spatial frequency do- 
main) and extracting the phase by using straightforward FFT coefficients 
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manipulations 



3.2 The L4Ci?£' phase-shift estimation: an introduction 

To introduce the I ACRE method ("Interferogram Analysis by Continuous 
Wavelet Transform Ridge Extraction"), let us observe that the sequence sz 
(and generally each sequence I{z = Z,x)) has the structure of a frequency- 
modulated sinusoid plus some corrections (noise, slowly varying background). 
It is therefore natural to try to extract the sz phase-shifts by using CWT 
techniques, with Ridge detection playing a relevant role. 

Consider the CWT map of the sequence sz@Z = 400 (see Fig. H). We 
can then try to apply the Ridge-extraction technique to the CWT map of sz 
to both denoise the sequence and extract the phase for each pixel position 
X. The Ridge sequence will be constituted by only the frequency-modulated 
components of sz, so that noise and background will be automatically dis- 
carded. This is the case for sz @Z = 400, as it is clear in Fig. ^. The phase 
sequence (pzix) for the analyzed array sz is then simply computed as the 
phase of the complex sequence of CWT at the Ridge: 

4>sz{x) = phase {(Wsz){x,aR{x))) , (6) 

and the phase-shift A0s^(x) is obtained as 

^<f>sz{x) = (f)sz{x) - M^) ^ (7) 

where (po{x) = kpX and kp is the wavenumber of the not-perturbed fringes. 

3.3 The I ACRE method step- by-step 

Let us now examine the recipe for the lACRE algorithm. Let I{z,x) be 
the gray-level image matrix of dimension M x N. The first steps are the 
estimation of the unperturbed fringe wavelength kp and (eventually) image 
filtering to slightly reduce noise. Next, for each Z G [1 M] we consider the 
sequence 

szix) = I{z = Z,x) ,x e [IN] 

and: 

• Compute the (complex) CWT map Wz{a,b) with the Morlet 
base in the Log sampling. To do this one must choose the number of 
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voices per octave A^^,. A large {N^ > 12) should be preferred if fast 
changes in the local frequency Q{x) are expected. In addition, if we ex- 
pect that in some regions the local frequency Qz{x) could have abrupt 
changes (local irregularities, structures, edges ...), a higher spatial res- 
olution is preferred {uq = 2tt, r < 1), while for regular behaviour (like 
the one of interferogram Int 1) a medium space-frequency resolutions 
should be used (r = 1). 

• Detect the (complex) Ridge sequence Rz{x) = Wsz{x, aji{x)) . 

• Compute the phase of Rz: 

4>z{x) = phase{Rz{x)) . 

• Estimate the phase-shift at z = Z as 

S(j){Z, x) = 4>z{x) — kp .X 

The result is a phase-shift matrix S(f){z, x) of dimension M x N. Phase 
unwrapping algorithms are then applied to the phase-shift map to eliminate 
unphysical phase jumps (this is the case for FFT-based results too). 

3.4 Comparison between the I ACRE and FFT-based 
performances 

We now apply the CWT-based and FFT-based methods to both real and 
simulated interferogram images. To start, we apply the lACRE method to 
the whole interferogram of Fig. |l|, which is corrupted from noise and shows 
strong reduction of fringe visibility and the presence of small scale periodical 
structures not related to the plasma properties. 

The lACRE output result is obtained with the image partially filtered 
from noise with a Median-Filter of mask size 3x3 pixels followed by a 
Wiener-Filter 5x5 and using A^^, = 12 voices per octave (low A^^,). The 
phase shift map is shown in Fig. ^ which should be compared with the 
FFT-method phase-shift of Fig. ^ obtained with the same filtered image. 
As it is clear from figures Figg. || and ||, while the CWT output seems to 
be accurate, FFT output is noisy and not free from unphysical phase jumps 
near the target, where a strong reduction of fringe visibility is present. 
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The higher accuracy of the MCi^E" method with respect to the FFT-based 
one is a very important characteristics of our new procedure. It enables an 
accurate search of small non- uniformity of the phase-shift map which are 
important to detect the growth of plasma instabilities as filamentation and 
self-focusing. 

To better check this point, we numerically build-up one interferogram in 
which we simulate the phase shift produced by a slowly-varying background 
plus some small scale filaments. Noise and reduction of fringe visibility are 
finally added to the interferometric image to better match the real interfer- 
ograms characteristics. 

The interferogram of Fig. |^ simulates a plasma with a background of 
maximum electronic density {fie/n^^Max = 0.1 with a Gaussian profile in the 
radial direction (with radius 75 /im) which is exponentially decreasing in the 
X direction. Three filaments are then added in different positions, each one 
with Gaussian density profile: 

= aexp(— +y)/r) 

ric 

with maximum density perturbation and radius [a = 0.005, r = 10 fim), [a = 
0.005, r = Sfim) and (a = 0.005, r = 6/im), respectively. Since the electronic 
density is everywhere much lower than the critical density, the linearity of the 
phase map with respect to the density is respected. We can then compute 
the perturbation of the phase-shift map in 2tt units {(f)2n = 0/(27r)) with 
respect to the background as 

1 r 5n 



6{A<f)2.) = -77^ / —dy 



whose maximum value is 



1 T 

S{A<p27T)Max = T^VTTa— , (8) 
Z An 



that is 5(A027r)Max = 0.18, 5(A027r)Max = 0.14 and S{A(j)2n)Max = 0.11, 

respectively. To detect these structures, the noise level of the phase-shift 
map should be a fraction of S{A(j)2n)Max- If cr(a^) is the standard deviation 
of the noise of each sequence of A02,r(-2, x) at x fixed, we could detect these 
structures if their amplitudes are for instance at "two sigma" with respect to 
the noise, that is cr(a;) < 0.09, cr^x) < 0.07 and cr^ < 0.05, respectively. The 
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standard deviation cr(x) (or one fraction of cr(x)) could be then be considered 
as a rough estimation of the "detectable threshold" in the phase-shift map". 

To estimate the accuracy of the phase-shift maps obtained with the 
I ACRE and FFT-based methods, we compute the phase-shift maps with 
the two methods (see Figg. ^ and ^ and we compare them with the known 
"true" phase map. We start the analysis by comparing some line-out of the 



two phase maps with the known simulate map. In Fig. |10] it is clear that 
the accuracy in the two phase methods is comparable in regions of the in- 
terferogram with low phase-shift, while for large phase-shifts the FFT-based 
output clearly fail in producing an accurate phase map. 

Denoting with A0^^ and A^^/'-^ the phase-shifts maps obtained with 
the two methods and with A0™™*^ the simulated phase map, we estimate 
the error map as the differences: 



Scwt{z,x) ^ A0^7^(z,x)-A0™™^( 



Z,X) 



£FFTiz,x) ^ A0^/^(;^,x)-A0™(;^,x), 

so that the sequences of the standard deviations of the noise in the phase 
map can be estimated as 

apFrix) = std{SFFT{z,x)) , 

c^cwt{x) = std{£cwTiz,x)) , (9) 

where std{f{z)) is the standard deviation of a sequence f{z). 

In Fig. is shown the behaviour of the error in both the I ACRE and 
FFT-based maps ((a)) while in b) the ratio R{x) between aFFT^x) and 
(Tcwt{.x) sequences is reported. The analysis of these figures confirms the 
claim that in small phase-shifts regions the MCi^E' method exhibits a slightly 
higher precision than the FFT-based one (the (Tfft / o^cwt sequence is about 
2), while in large phase-shift regions the sensibility of the I ACRE method is 
much higher than the one of the FFT-based one. For example, assuming the 
sequence a{x) as an estimation of the phase-shift sensibihty, since for gaus- 
sian density profiles (5(A(/)2^)Max is proportional to the structure radius r and 
the maximum density perturbation a (see H), the sequence R = (Tfft I c^cwt 
could be interpreted as a rough estimation of the ratio between the mini- 
mum product ar detectable with the FFT-based and MCi^E" techniques: 

ocwT ^ (ar))^^^"^- 
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Futhermore, since for the I ACRE method the (Tcwt sequence is every- 
where below the value 0.03 (see the "detection thresholds" reported above), 
we are confident that all the three filaments could be detected. This is not the 
case for the FFT-based method output because in the region X G [220, 250] 
the apFT sequence is in the range 0.05 — 0.15, which is over the minimum of 
the detection thresholds. 

We conclude the analysis of the interferogram of Fig. 7 by checking the 
behaviour of the phase maps when an algorithm for the automatic extraction 
of small scale perturbations is applied to the phase-shift maps. The algorithm 
utilized is very simple and consists of two main steps: 

• The decomposition of the map /S.(j)2-K in a 'Large scale' component (the 
background) ^<t>2-K ^^nd a 'Small scale' component 5(A027r) (the struc- 
tures + noise) by using a Smoothing B-spline fitting for each line-out 
of the phase map at x fixed. 

• The filtering of the small scale component 5(A</)27r) with a "two sigma" 
cutoff. As explained before, provided that structures in the (5(A027r) 
give a negligible contribution in the Root-Mean-Square of the map, the 
standard deviations apFrix) and (Tfft{x) can be computed as 

acwAx) = std{5{M'i^^){z,x)), (10) 



In Fig. [T2| the filtered at "two sigma" 'Small scale' phase maps obtained 
with the two methods are reported. As expected, the filtered map of the 
I ACRE method clearly shows the presence of the three filaments, while in the 
FFT-based map some regions of the map with strong presence of noise could 
be interpreted as false small scale structures so no clear filaments detection 
is possible. 



4 Conclusions 

With the help of one real and one simulated interferograms we showed that 
the I ACRE method is more accurate and robust than the FFT-based one. For 
the simulated interferogram the smallest detectable phase-shift perturbation 
(with respect to the background) obtained with the I ACRE method is in 
the mean 0.5 times the one obtained with the FFT-based one, with possible 
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further decrease in higher density regions. In addition the outputs of I ACRE 
arc free from unphysical phase jumps both in the real and the simulated 
interferograms, while FFT-based map is in both cases affected by a large 
region near the target where phase jumps cannot be removed by conventional 
unwrapping procedures. The higher robustness and sensibility of the I ACRE 
method can be addressed both to the wide adaptabihty of the CWT tool 
to the actual image and the intrinsic strong noise suppression in the Ridge 
Extraction procedure. 
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Figures Caption 



Fig. 1 A sample interferogram of a plasma produced by laser explosion of a 
0.5/im thick, 400/im diameter Aluminium dot coated onto a 0.1/im plastic 
stripe support. The interferogram was taken, perpendicularly to the strip 
surface, 3.0ns after the peak of the plasma forming pulses using a modified 
Nomarski interferometer. The intensity on target was 8.5 x lO^^W / 071"^ . The 
probe pulse-length was lOOps and the probe wavelength was 0.53fim. For 
details on the experimental set-up see 0. 

Fig. 2 Line-out of the fringe intensity (sequence sz) of the interferogram 
of Fig. 1 at Z = 400. 

Fig. 3 The real part and absolute value of the CWT maps of the signal 
sz@ Z = 400 (interferogram of Fig. 1). 

Fig. 4 The sequence of real part of the Ridge sequence of signal sz@Z = 
400 (interferogram of Fig.l). 

Fig. 5 The phase-shift map (in 2tt units) obtained from the interferogram 
of Fig. 1 after suitable filtering. lACRE method. 

Fig. 6 The phase-shift map (in 27r units) obtained from the interferogram 
of Fig. 1 after suitable filtering. FFT-based method. 

Fig. 7 A simulated interferogram of a plasma containing three small 
filaments. 

Fig. 8 Phase-shift map of the simulated interferogram of Fig. 7. lACRE 
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method. 



Fig. 9 Phase-shift map of the simulated interferogram of Fig. 7. FFT- 
based method . 



Fig. 10 Line-out of the phase-shift maps of the simulated interferogram 
of Fig. 7. The lACRE and FFT-based methods outputs are compared 
with the 'true' simulated map. 

Fig. 11 (a) Standard deviations of the error in the phase-shift map of 
interferogram of Fig. 7 computed via lACRE and FFT-based methods, (b) 
Ratio between the standard deviations of the error in the phase-shift map of 

interferogram of Fig. 7 computed via FFT-based and lACRE methods. For 
density perturbations 5n with gaussian density profile in the z direction of 
amplitude a and radius r, the sequence (ypFT /(^cwt can also be interpreted 
as the ratio between the minimum product ar detectable with the FFT- 
based and lACRE methods: (Tfft/(7cwt ~ {ar)fFT/{c(r)^WT- 

Fig. 12 Filtered map at "two sigma" of the 'Small scale' component of 
the phase-shift map of interferogram of Fig. 7. (a) FFT-based method: two 
filaments could be detected but other unreal structures survive to the "two 
sigma" filter, (b) //ICi^E method: three filaments are clearly detected. 
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